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Abstract 



We propose a formulation of the long-distance dynamics of gauge theories 

at finite temperature on a lattice in Minkowski space, including the effects 

of hard thermal loops on the dynamics of the long wavelength modes. Our 

approach is based on the dual classical limits of quantum fields as waves 

and particles in the infrared and ultraviolet limits, respectively. It exhibits 

manifest invariance under space-dependent lattice gauge transformations and 

conserves Gauss' law. 

02.60.Cb, ll.15.Ha, ll.15.Kc, 12.38.Mh 



Typeset using REVTgX 



*chaoran@phy. duke.edu 
^muller@phy.duke.edu 



1 



In the past few years great efforts have been made in perturbation theory to calculate 
transport properties of thermal gauge fields from the low-energy effective action including 
the contribution of hard thermal loops Unfortunately, perturbative calculations fail 

in some important cases such as the damping of a traveling mode in the QCD plasma ||, 
color conductivity ]7j|, and winding number diffusion |3j due to the presence of singularities 
associated with the static magnetic gauge sector. This difficulty has motivated numerical 
simulations of the dynamics of classical gauge fields in Minkowski space [|9Hl2|l, which were 



based on the Hamiltonian formulation of gauge theory on a spatial lattice [13[ . 



These studies have been criticized fll^Jlql , because they do not properly account for the 
influence of hard thermal loops which modifies the dynamics of the long-distance modes 
of the gauge field. For scalar field theories there exists a straightforward remedy ||i4| , |T6 



One introduces a momentum cut-off k c and describes the influence of the high-momentum 
modes, which are essentially quantum mechanical, in perturbation theory. The dynamics of 
the long-distance modes then becomes dissipative and noisy fl6| , [T7|j . In the case of gauge 
theories, gauge invariance dictates the use of a lattice discretization eliminating modes with 
eigenvalues of the kinetic momentum operator (— zV — gA) larger than it /a, where a is the 
lattice spacing. However, the construction of the appropriate lattice action for the soft fields 
involving hard thermal loops encounters technical difficulties [|14| . 



Recently, Huet and Son |L8| have argued that it is possible to construct an effective 



classical dynamics for quasi-static, long wavelength modes of the gauge field containing 
dissipative and stochastic terms. The resulting equations are valid in the limit ui <C k < g 2 T 
where uo is the frequency and k the wave vector. Even in this limit the noise term is strongly 
nonlocal and no efficient numerical treatment was proposed by the authors. In the general 
case, for u,k <C T, the hard thermal loop action contains an infinite number of spatially 
and temporally nonlocal, dissipative and stochastic vertices, which are difficult to treat 
numerically. 

We propose to circumvent these difficulties by representing the hard thermal modes as 
classical colored particles propagating in the background of the soft gauge fields. Explicitly 



treating these modes in terms of classical particles, rather than integrating them out, leads 
to a set of local dynamical equations which can be efficiently solved by numerical intergration 
after lattice discretization. Below we will construct a lattice version of these equations that 
is gauge invariant and conserves Gauss' law. Because the dynamics of the hard modes is 
treated explicitly in this formulation, it is not necessary to assume that they are thermally 
populated. One well, consider dynamical situations far off equilibrium, where the 

density matrix of hard modes of the gauge field is characterized by some large scale. Such 
conditions are, for example, of interest in the context of equilibration processes occurring in 
high energy nuclear collisions. 

The representation of the high momentum components of the gauge field in terms of 
classical particles propagating in classical gauge fields can be justified by the eikonal limit 
of the Yang-Mills equations. Heinz |19j and Kelly et al. [2(J have shown that a thermal 



ensemble of particles obeying Wong's equations |2T| generates the correct hard thermal loop 



action |T||J] for soft gauge fields. The proof is based on linear response theory, combined 
with the explicit gauge covariance of the classical equations. In the following we briefly 
review the continuum formulation of the classical transport theory and then show how 
the equations can be implemented on a spatial lattice. 

The classical transport theory for nonabelian gauge fields starts from the Boltzmann 
equation 



if 



" 9QaF ^dp~u ~~ gUhcA ^ QC dQ a 
together with the Yang-Mills equations 



f(x,p,Q) = C[f] (1) 



D,F^ = g J [dpdQ]p v Qf{x,p, Q) = f(x) , (2) 

where f(x,p, Q) denotes the one-particle phase space distribution of classical particles, Q is 
the classical nonabelian charge carried by the hard gluons, and is the gauge covariant 
derivative. Note that p represents the kinetic (not the canonical) momentum of the particles 
and therefore is gauge invariant. 



For our purposes, it is sufficient to consider the Vlasov limit, neglecting the collision 
term C[f). The transport equation ([!]) is then solved by an ensemble of test particles 

f{x,p,Q) = ±r y E, s (x-£i)S(p-Pi)HQ-Qi) > ( 3 ) 

where N is the total number of particles. The space, momentum, and color coordinates of 
the particles obey Wong's equations 

£i = Vi, (4) 
p i =9Qf[E a (^+v i xB a {^] , (5) 
Q i = -igvnA^ i ), Qi] . (6) 

The index i enumerates the particles and £•* = (t, = (1, v) with Vi = Pi/\Pi\ being 

the velocity of the ith particle. Note that i> M , as defined here, is not a Lorentz four- vector. 
A dot denotes a time derivative. denotes the vector potential, and E,B are the color 
electric and magnetic fields, respectively. The right-hand side of ([5]) is a generalization of 
the electrodynamical Lorentz force. Wong's equations can be written in manifestly covari- 
ant form, but here we have chosen a representation that is convenient for the numerical 
implementation in the context of a Hamiltonian lattice gauge theory. 

The Hamiltonian equations of motion for the gauge field and Gauss' law are obtained 
from (HD as the components v = 1,2,3 and v = 0, respectively. The charge current j M = (p, j) 
is defined as 

f(x)=gY.Qr<5\x-£ l ) ■ (7) 

i 

Since E, B, and Q transform covariantly under gauge transformations, it is easy to 
see that (f|) and (|5|) are gauge invariant. The gauge covariance property of (||) is best 
demonstrated by first recognizing that (||) has the following formal solution (omitting the 
index i): 

Q(t) =U(t,0)Q(0)U\t,0) , (8) 
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where U(t, 0) is the parallel transport operator along the particle's world line: 



U(t,0) = Texp 



(9) 



U obeys the equation 



^^(OW(t,0) = 0, W(0,0) = 1. (10) 



Under a gauge transformation G, 14 (t) transforms as U(t,0) — > G(x(t))U(t,Q)G^(x(0)). 
Hence (U) is gauge covariant, and so is (§). It is worth noting that both (H) and (§) conserve 
the magnitude of the charges: 

As the charged particle travels in the background field, its charge rotates in color space with 
an angular velocity u a = gv^A^^), consistent with the notion of the gauge field A 1 as the 
connection on a curved manifold on which the charge of a moving particle undergoes parallel 
transport. 

We now describe how the equations (|2],f|,|5|||) can be solved numerically. We discretize (g) 
on a lattice with lattice spacing a. Then the soft gauge fields are restricted to modes with 
k < k c = Tc/a. To avoid double-counting of modes the hard particles will be restricted to 
kinetic momenta p > k c . The precise connection between k c and a can be established by the 
requirement that the plasmon mass takes on its correct value found in the continuum theory: 

= y/T^gT '/3. In the proposed simulation, we are mainly interested in modeling physics at 
the energy scale of g 2 T and the dominant scattering process between hard particles happens 
at the scale of gT. In order to get these right on the lattice, we must set up the simulation 
in such a way that the following relation between the length scales is satisfied: 

a < (gT)- 1 < (g 2 T)- 1 < Na , (12) 

where Na is the lattice size. 

For discretization of the fields, we use the Kogut-Susskind (KS) model, which is summa- 
rized below. Improved lattice Hamiltonians are also possible |[23|| . In the KS scheme, one 



chooses the temporal gauge A = and expresses the gauge field in terms of variables U Xji 
associated with the link (x, i) directed from a site x to its nearest neighbor x + i: 

U Xti = exp [-igaAt(x) r a /2] = U\ +i ^ . (13) 

Under a gauge transformation G, the link variables transform as 

U x ,^G{x)U x ,G\x + i) . (14) 

A plaquette variable is defined as the product of four link variables associated with the sides 
of an elementary plaquette (x,ij): 

U Xj ij U x ^U x j r j i ^jU x j r i^.j^iU x j r j^j . (15) 

The links are directed and hence the plaquettes are oriented. The electric and the magnetic 
fields are defined as 

E x ,i = —U x ,iU x ,i , (16) 
iga 

Bx,k = 4^q2 tki 3^x,ij — U x> ij) , (17) 

where E x>i is associated with the link (x, i) and B x k with the plaquette (x,ij). Both E x>i 
and B X) k transform covariantly under a gauge transformation G(x). 

In the spirit of the Hamiltonian formalism, we choose U Xji and E Xji as the basic dynamic 
variables. They obey the following equations of motion: 

U Xti = igaE Xji U X:i , (18) 

^,< = ^E(<«-^-)-j.,<, (19) 

and are subject to the constraint of Gauss' law: 

- E \E Xii - Ul_ hl E x . hl U x .^ - p x = . (20) 
a i 

In order to define the charge current four-vector j£ on the lattice, we take each site x 
as the center of a cubic cell C x of size a 3 . The color charge of every particle in C x will be 
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counted as contribution to the charge density p x . Any particles entering or leaving the box 
during a given time step At will contribute to the component of the color current normal to 
the face of the cube that is penetrated p^j : 



p»(f) = 4 E Qk® , (2i) 



a 3 



(i) 

j x4 (t + At/2) = -J- £ [Q k (t) -Q k (t + At)] , (22) 
a 1X1 kaC x 



where the notation in ( p2|) indicates that only particles entering or leaving along the link 
connecting x and x + i are counted. 

The variables to be advanced in time according to their equations of motion are U x j, 
E x ,i, Pki an d Qk- This can be achieved with a modified leapfrog algorithm that conserves 
energy and satisfies Gauss' law. It is most convenient to choose an update scheme in which 
the link variables U x ^ are defined at half integer time steps while E Xji , £ fc , p k , and Q k are 
defined at integer time steps [p5|1 . We only briefly discuss the momentum update here. 



The momentum gets contributions from both the electric and the magnetic fields. In the 
same spirit as our definition of the current j x j in (p2]), only those particles which move from 
one cell to another during one time step obtain a momentum kick by the electric field: 



p Ki {t + At) = Pfc>i (t) + gQ a k {t)E a x M) 



1 tr [Q k (t)Q k (t)} ' 
v kii { 2a 2 tr[Q k (t)E X:i (t)] 



(23) 



if k G C x at time t and k G C x+ i at time t + At. The expression in the brackets can be 
regarded as the effective time needed for a particle to transit from one cell to its nearest 
neighbor. This choice balances the change in the energy of the particle and the change in the 
energy of the lattice fields incurred by the transition of that particle from one cell to another. 
The influence of the magnetic field on the particle momenta is updated continuously during 
each time step. Such a modified leapfrog algorithm is stable for At/a ~ 0.1. 

Finally, we point out that the definition (p2] ) ensures the conservation of the color current 
during each time step: 
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^ (p x (t + At) - Px (t)) + ^Y1 (*M ~ '', ,,./-• ,,1'r ,,) = , (24) 

i 

where the quantities in the sum are defined at t + At/2. This, together with the update 
scheme described above, automatically ensures the exact conservation of the left-hand side of 



Gauss' law (|20D from one time step to the next. An important feature of our implementation 
is that all equations transform correctly under spatial gauge transformations. 

The dynamical evolution of (soft) gauge field and (hard) particles described by the equa- 
tions of motion is essentially classical. However, we need to incorporate certain quantum 
features into the simulation, which are embodied in the initial conditions. There are two 
places where quantum physics enters. 

First, while the energy and the momentum of a charged particle are continuous classical 
variables, the nonabelian charge Q a contains a factor of h and has a fixed magnitude. In 
the classical limit the color charge of a gauge boson rotates in color space like a three-vector 
with fixed length just as the spin of an electron rotates in a magnetic field. In the quantum 
case, the charges Q a are g-numbers and obey the SU(2) Lie algebra: 

[Q a , Q h ] = ihe abc Q c . (25) 

In the semiclassical limit, we treat them as c- numbers but retain their magnitude as pro- 



portional to the Casimir operator in the adjoint representation [19 



]T Q a Q a = 2h 2 , (26) 

a 

which is conserved by (^|). For particles in the fundamental representation of SU(2), such as 
fermions, the right-hand side would be replaced by 3h 2 /4. 

Second, while the dynamics is classical, we require that both the particles representing 
hard thermal gluons and the fields obey Bose statistics. To illustrate this point, we consider 
the initialization of the particle ensemble at a certain temperature T. According to Bose 
statistics, the particles should be initialized with the distribution 



f(x,p,Q) = 6 ^ {p) n ^ 5(Q 2 -2h 2 ), (27) 



where e(p) is the energy of a particle with momentum p. The initial particle number density 
is then given by 

^^^/t^i^-^' (28) 

where the factor 2N C counts the spin and color degeneracies. The linear response of the 
distribution f(x, p,Q) to the classical field then reproduces the HTL polarization function 

MM- 
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